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Phase transitions in a gas of anyons 
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We continue our numerical Monte Carlo simulation of a gas of closed loops on a 3 dimensional 
lattice, however now in the presence of a topological term added to the action corresponding to 
the total linking number between the loops. We compute the linking number using certain notions 
from knot theory. Adding the topological term converts the particles into anyons. Using the corre- 
spondence that the model is an effective theory that describes the 2+1-dimensional Abelian Higgs 
model in the asymptotic strong coupling regime, the topological linking number simply corresponds 
to the addition to the action of the Chern-Simons term. We find the following new results. The 
system continues to exhibit a phase transition as a function of the anyon mass as it becomes small 
[T], although the phases do not change the manifestation of the symmetry. The Chern-Simons term 
has no effect on the Wilson loop, but it does affect the 't Hooft loop. For a given configuration it 
adds the linking number of the 't Hooft loop with all of the dynamical vortex loops to the action. 
We find that both the Wilson loop and the 't Hooft loop exhibit a perimeter law even though there 
are no massless particles in the theory, which is unexpected. 

PACS numbers: 11. 15. Ha, 11.15.-q, 11. 15. Ex, 04.60. Nc, 02.70.Ss, 05.30.Pr 



I. INTRODUCTION 

In this paper, we continue our study of a 3-dimensional 
lattice loop gas [IJ, adding a topological interaction which 
counts the total linking number of the loops with each 
other. This Euclidean theory corresponds to a 2+1- 
dimensional model where the particles and anti-particles 
are non-interacting hard-sphere anyons. Thus the only 
interactions they suffer are a statistical, Aharonov-Bohm 
[3j type interaction and an infinite short range repulsion 
between particles (and anti-particles) which does not al- 
low them to come close to each other. The Aharonov- 
Bohm interaction induces a change of phase of the wave 
function when one particle encircles another, giving a 
phase dependent on the coefficient of the topological 
term. Indeed, it is easy to picture a process whereby two 
particle anti-particle pairs are created, live for a while 
and then annihilate. If the two particles make a 360° 
turn about one another, the process is described by two 
loops which link once with each other. The topological 
term contributes a phase e*" to the corresponding ampli- 
tude, where k is the coefficient of the topological term in 
the action. 

The theory admits an interpretation as an asymptotic, 
strong coupling, effective Euclidean lattice description of 
the 2+1 dimensional Abelian Higgs model [1_ in the sym- 
metry broken sector, the topological interaction corre- 
sponding to the addition of a Chern-Simons term (CS), 
with coefficient k/2. We will lean heavily on our pre- 
vious work pL,; thus we refer the reader to that article 
for our conventions. In the symmetry-broken phase, in 
the asymptotically strong-coupling limit, the theory de- 
scribes the dynamics of non-interacting vortices and anti- 
vortices which now behave as non-interacting anyons due 
to the Chern-Simons term. The system continues to ex- 



hibit the phase transition that was studied in [Tl. 

The Euclidean theory has an action that is no longer 
real, since the Chern-Simons term adds an imaginary 
term to the action. This is in principle an impediment 
to the simulation of the theory using the Monte Carlo 
method since the exponential of the Euclidean action can- 
not be interpreted as a probability density. However, we 
can continue to use the real part of the action to give us 
the probability distribution. It has the same symmetries 
as the full action, although it is additionally symmet- 
ric under parity. Then the Chern-Simons term is simply 
a bounded unimodular phase which can be integrated 
against the measure that is defined by the real part of 
the action. 

The Chern-Simons term for the configurations that we 
are left with (non-intersecting closed loops on the lattice) 
is simply equal to twice the total linking number of all 
the loops [2J. We can see this quite easily by observing 
that since the Chern-Simons integrand is proportional to 
the magnetic field, the integral reduces to a set of line 
integrals along the vortex fiux lines: 
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In the second line the integral can be separated into a 
sum of integrals along each vortex line, denoted Ci, and 
a two-dimensional integral in the transverse direction. 
The transverse integral is simply equal to the total flux 
in the vortex line, which is 27r, while the integral of the 



gauge field along a given vortex line just gives, via Am- 
pere's law, the total flux of all the other vortex lines that 
link with the first vortex line. Thus each integral in the 
sum is given by 27riVi(Ci), where NlICi) is exactly the 
linking number of vortex line Ci with all of the other 
vortex lines. Summing over all the curves clearly gives 
twice the total linking number of all the vortex loops, 
TliC ^L{Ci) = 2Nt; where Nt is the total linking num- 
ber of the configuration of vortex loops; i.e.. 
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It is this total linking number that we must compute. 



II. COMPUTING THE LINKING NUMBER 

The closed non-intersecting loops were generated on a 
body centered cubic (bcc) lattice of size 100'^ by plac- 
ing the cube roots of unity randomly on the vertices, 
with all points on the surface of the lattice assigned the 
same value. The lattice can be thought of as filling space 
with (non-regular) tetrahedra. Each cube contains six 
pyramids; adding the diagonal of the cubic sides in a sys- 
tematic way throughout the lattice divides each pyramid 
into two non-regular tetrahedra (see pi for details). If 
the change in phase of the cube roots of unity around 
one of the triangular faces of a tetrahedron is equal to 
±27r, we say a length of vortex fiux tube has entered or 
exited the tetrahedron through that face. Suppose the 
fiux entered the tetrahedron. It is a quick exercise to 
conclude that whatever cube root of unity is placed on 
the fourth vertex of the tetrahedron, the fiux must exit 
the tetrahedron through one of the other faces. However 
then it enters another tetrahedron since the tetrahedra 
fill space, and the analysis must be repeated. The loop 
must close, since the boundary condition used means no 
surface triangle has a fiux passing through it. The loops 
so defined exist in the dual lattice to the inital tetrahedral 
lattice. 

To compute the linking number of any given loop with 
all of the other loops, we need to simply compute the 
fiux that passes through the given loop, since each other 
loop that links with it carries one unit of fiux. To com- 
pute this fiux is actually quite difficult since the given 
loop and in fact all of the loops exist in the dual lattice 
where none of the original variables are defined. On the 
other hand, the fiux that passes through a loop that is 
defined along the links and vertices of the original lattice 
is trivially calculated: we simply calculate the change of 
phase of the cube roots of unity as we pass through the 
vertices of such a loop. Thus if we can systematically de- 
form the given vortex loop on the dual lattice to a loop 
on the original lattice, the calculation is straightforward. 
This deformation is a virtual deformation through space: 
we do not alter the configuration of the cube roots of 
unity. Hence, the calculation of the linking number via 



the fiux that passes through the deformed loop also adds 
the linking number of the deformed loop with the orig- 
inal loop, what is reasonably interpreted as a definition 
for the self-linking number of the original vortex loop. 
We must subtract off this self-linking number from the 
calculation of the linking number of the deformed loop in 
order to obtain the linking number of the original loop. 
The deformed loop and the original loop form a ribbon, 
and it is the linking number of the two curves that form 
the edges of the ribbon that we must calculate. This 
is normally called the self-linking number of the original 
loop. What we need to do is to subtract off this self- 
linking number. The self-linking number of a loop is not 
a completely well defined quantity |1] : it depends on the 
deformation that gives rise to the ribbon. However in 
this case the deformation is defined by requiring that the 
deformed loop live on the original lattice. The self-linking 
number satisfies the relation [Q 



Self linking number = Twist -I- Writhe. 



(3) 



The Twist of the ribbon is intuitively defined as the num- 
ber of times the orthonormal frame, defined by tangent 
vector, the perpendicular displacement vector in the in- 
stantaneous plane of the ribbon and the mutually orthog- 
onal vector product of the preceding two vectors, twists 
around the original curve as you go around the loop. Sur- 
prisingly, the Twist is not an integer. The Writhe is a 
more complicated object which also is not an integer, see 
[6] for details. It is intuitively the number of coils with- 
out twist that the ribbon suffers. Of course the sum of 
the two is an integer, the self-linking number of the rib- 
bon. A familiar example might correspond to a coiled, 
untwisted extension cord. When we straighten it out, in- 
variably it becomes all twisted up; this is simply because 
writhe is converted into twist, as the self- linking number 
is essentially invariant. 

However, the self-linking number is in fact the ordinary 
linking number of the deformed loop with the original 
loop. The Gauss formula for the linking number is given 
as two line integrals over the locations of the two loops 
which link: 
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where dx and dy are the line elements along the two 
loops respectively. These integrals are prohibitively time 
consuming to calculate directly since the loops are polyg- 
onal with of the order of 10^ segments. Hence one is left 
with about 10^^ integrals to do, which, although not im- 
possible, is not feasible with our computing resources. 
However there is a simpler way to compute this linking 
number, using knot theory [5_. 

The idea is the following. We can project the knot 
onto a two dimensional plane, keeping track only of the 
the sense of the crossings of the segments of one loop 
with the other in the projection. Then a simple sum 
of the association of ±1 to the crossings, depending on 



which segment of which loop is on top of the other and 
its direction, gives the linking number. The important 
point is to choose the direction of projection that will 
yield the most simplifying two dimensional projection. 
Indeed, an arbitrary projection will yield just as many 
independent crossings as there are links in the two loops, 
not significantly reducing the number of calculations (see 
Fig. §. 




FIG. 1. (color online) An example of a vortex loop (thin, 
yellow) and its deformation to the lattice (thick, red) viewed 
from a random direction. 



However, projecting along the diagonal of the lattice 
(the (1,1,1) direction in coordinate space) actually yields 
exactly a regular triangular lattice on the projected two 
plane. The original loop, which passes through the dual 
lattice, projects to the dual lattice of the two dimensional 
triangular lattice, while the deformed loop, of course, 
projects directly to the links of the triangular lattice (see 
Fig[2f. 

Thus the crossings are unambiguous and occur at a 
small, finite number of intersection points. It is easy to 
keep track of the segments of each loop, and their rela- 
tive heights. This simply amounts to a re-indexing of the 
data which is already stored in the computer, in a new 
system of coordinates given by the triangular lattice in 
the two dimensional projected plane and the height along 
the (1, 1, 1) direction. The original calculation of the link- 
ing number which would take several days of computer 
time reduces to a few seconds. 




FIG. 2. (color online) The same vortex loop and its defor- 
mation as in Fig. [Ilviewed along the (1, 1, 1) direction. 



III. CALCULATING IN THE CHERN-SIMONS 
THEORY 

The Monte Carlo process for generating a set of con- 
figurations in the presence of the Chern-Simons term is 
not straightforward. The point is that the Chern-Simons 
term adds an imaginary contribution to the Euclidean 
action; hence it cannot be used to define the probabil- 
ity distribution required to obtain the equilibrium con- 
figurations via the Monte Carlo method. However, we 
obtain the set of equilibrium configurations by using the 
Boltzmann weight given by the total length of the vortex 
loops. Thus the Boltzmann weight is unaffected by the 
Chern-Simons term. 

This is a classic problem for numerical simulations; it 
arises in theories with fermions where it is called the 
sign problem, but also in the context of topological terms 
which are odd under time reversal; a fuller explanantion 
of this is provided, for example, in [TU]. The solution in 
the present case is based on the following logic. The (Eu- 
clidean) Feynman path integral instructs us to integrate 
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over the space of field configurations, with appropriate 
boundary conditions, to calculate a quantum amplitude 
(O). However, the actual measure on this space is not 
specified. There does not exist a naive, canonical, trans- 
lationally invariant (translations in the space of func- 
tions) measure on the space of functions [TT] that we 
can use. However perfectly well defined Gaussian mea- 
sures do exist, and are defined by quadratic (free) field 
theories. It is the problem of constructive quantum field 



theory to prove that the set of ampHtudes created by the 
Feynnian path integral using the Gaussian measure on 
function space, integrating the exponential of the inter- 
actions in the Lagrangian against this measure, provides 
a full set of finite Greens functions that satisfy the Wight- 
man axioms |12| , which then allow for a reconstruction of 
the corresponding, interacting quantum field theory. In 
the present case, the interaction could indeed comprise 
of the Chern-Simons term, which we note, remains imag- 
inary in Euclidean space. The Gaussian measure suffices 
to define the quantum field theory, modulo some amount 
of (infinite) renormalisations of coupling constants, how- 
ever, the net effect of the construction of the quantum 
field theory is in fact to define a measure corresponding 
to the exponential of the full action, not just the Gaus- 
sian part. There are several examples of this kind of 
construction of non-trivial quantum field theory |13| and 
the corresponding measures. 

The important property that the measure, used to con- 
struct the quantum field theory, must satisfy is that it 
respect all the symmetries of the theory that one is try- 
ing to construct. In our case, we do not use only the 
Gaussian part of the action, since we are doing a non- 
perturbative, lattice based numerical analysis. We use 
the full real part of the action with the approximation 
that yields, for the action, the total length of the vortex 
loops for the field configurations that we consider. The 
full real part of the action actually has more symmetry 
than the theory with the Chern-Simons term added since 
parity and time reversal are not preserved. However this 
does not cause any problem since all amplitudes are cal- 
culated with the Chern-Simons term inserted. 

Hence we use the full real part of the action to define 
the measure on the space of field configurations, and in 
the context of the Monte-Carlo method, the set of equi- 
librium configurations, say a total number N. After this, 
the Chern-Simons term simply gives a uni-modular phase 
that can be integrated against this measure. When cal- 
culating actual matrix elements of an operator, we must 
calculate the average of the operator with the phase com- 
ing from the Chern-Simons term inserted and then divide 
by the partition function, again defined with the same 
phase inserted. Explicitly we get 
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where Ci stands for the ith configuration. Normally, both 
numerator and denominator contain a factor of 1/iV; 
however, here it cancels between them. 



III.l. Chern-Simons term 

The Chern-Simons term in the functional integral gives 
exactly the linking number of all the dynamical loops. 
We use knot theoretic techniques to compute the average 
of the Chern-Simons term in the set of equilibrium con- 
figurations, for different values of the mass fj,. We see that 



the expectation value drops to zero remarkably quickly 
as the system passes through the transition at about 
/i — 0.152 for decreasing /i. The graph for the expec- 
tation value of the Chern-Simons term is given in Fig. [3] 
for various values of /i. This average of the Chern-Simons 



1,0 


■ r 

I.I. 


ill** **»»».^ ®®*®<f 




0,8 


Lattice size 102^ 

H 
a 0,13 


0,6 


- 


■ »* ♦♦ *j. 


® 0,135 
A 0,14 
+ 0,144 






■*.* ♦!* *♦ 


•J. 0,146 


0,4 


" 




4 0,148 
* 0,15 
e 0,152 


n? 




■*e * **:* 


* 0,154 






■■1- ♦* 
■■«•• 


« 0,158 




1 


c 0,16 




1 



FIG. 3. (color online) The average value of the Chern-Simons 
term as a function of k. 

term serves as the partition function when computing the 
expectation value of any operator. Clearly proceeding to 
values of K > 0.08 is not possible for /i < 0.15. 



III. 2. Wilson loop 

The Wilson loop [7] is defined as the expectation value 
of the operator 



ly = g-i(g/e)f Af.dx'^ 



(7) 



where e is the fundamental charge in the model and the 
integral in the exponent goes along a closed, fixed rect- 
angular path of width L and length T. For a given con- 
figuration Ci, the exponent in the Wilson loop is equal to 
the linking number of the curve defining the Wilson loop 
and the dynamical vortex loops in the configuration, say 
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The calculation of the expectation value of the Wilson 
loop can be done by first calculating the average value of 
the Wilson loop for fixed total linking number, and then 
performing the sum over these average values, weighted 
by the number of configurations with the fixed total link- 
ing number. Using the notation M^Nt) for the number 
of configurations with fixed total linking number Nt, and 
Ci^Nr ^s an index for these configurations, we have: 
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The term in parentheses in the numerator is the av- 
erage value of the Wilson loop with fixed total linking 
number. If this is independent of the value of the total 
linking number, then it comes out of the sum, and in 
fact the sums in the numerator and denominator cancel, 
yielding 
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which is in fact independent of k. This is exactly what 
we find with our numerical simulation. In Fig. [4] we plot 
the value of the Wilson loop for different values of /i as a 
function of k, at a fixed value lnqje = O.IStt. Evidently, 
the average of the Wilson loop does not depend on k for 
any value of p. for k < 0.08. 

There is an apparent dependence in the graphs for 
small 11, as K exceeds the value ^ 0.08. However, at 
this point, as we can see from Fig. [3] the average value 
of the Chern-Simons term, by which we divide, becomes 
very small, and we no longer trust the numerical results. 
Indeed, for small /x, the number of configurations at the 
extremities of the distribution of the total linking num- 
ber (i.e. at large total linking number, which only occurs 
for small /i) becomes only a handful. Then the average 
of the Wilson loop operator for these configurations de- 
viates wildly with respect to the average when there are 
many configurations (i.e. about 100) deviations which 
are magnified when the denominator also becomes small. 
We have verified that if we increase the number of con- 
figurations that we have available at fixed total linking 
number, then the average of the Wilson loop for this set 
of configurations converges to the k independent value as 
the number of configurations becomes large ^ 100. 



III.3. 't Hooft loop 

The 't Hooft loop [8_ corresponds to the insertion of 
a singular magnetic fiux tube along a contour of a fixed 
rectangular loop of width L and length T. It is the dual 
object to the Wilson loop P] . It is important to note that 
this is not a vortex loop, but just a gauge field loop. Thus 
the functional integral over the gauge fields is subject to 
the constraint that such a magnetic fiux loop exists at the 
given fixed position. The Monte Carlo method of gener- 
ating the equilibrium configurations is unchanged, using 
as before only the real part of the full action, with the 
(infinite) action of the 't Hooft loop subtracted off and 
with our strong coupling approximation. Then the equi- 
librium configurations are comprised of configurations of 
closed vortex loops appended by the 't Hooft loop. In the 
presence of the Chern-Simons term, the 't Hooft loop sim- 
ply adds kNuhl to the action, where Nuhl is the linking 
number of the 't Hooft loop with all the dynamical vortex 
loops. Hence the 't Hooft loop is given by the average 



CtH) 



Ec 



^iKNiffj -Se+iuNt 



y^ q-Se+ii^Nt 



(11) 



In Fig. [5] we plot the average value of the 't Hooft loop 
as a function of k, for various values of the mass /i. The 
points in the graphs beyond k = 0.08 should not be 
trusted for small values of p, since the errors are not 
under control. The average value of the Chern-Simons 
term in the denominator, as we apply Eqn. (|6|, becomes 
vanishingly small. We note that in contrast with the 
Wilson loop, the 't Hooft loop has a clear dependence on 
the coefficient of the Chern-Simons term. The 't Hooft 
loop is constant (equal to 1) in the absence of the Chern- 
Simons term, but is a function of the coefficient of the 
Chern-Simons term, in its presence. 
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FIG. 4. (color online) The Wilson loop for various values of 
the coefficient of the Chern-Simons term: k. 



FIG. 5. (color online) The average value of the 't Hooft loop 
in the presence of a Chern-Simons term as a function of n. 



IV. DISCUSSION AND CONCLUSIONS 

The main result that we find is that the Wilson loop 
[7] is independent of the Chern-Simons term while the 
't Hooft loop [8] is not. Both are expressed as the ex- 
ponential of a contribution to the free energy. In fact 
for the 't Hooft loop the contribution to the free energy 
vanishes in the absence of the Chern-Simons term but is 
proportional to its coefficient in its presence. The free 
energy for the Wilson loop is simply independent of the 
Chern-Simons term. 

These results are already surprising; however, there is 
an even more astonishing result. The fact is that both 
the Wilson loop and the 't Hooft loop have perimeter 
law behaviour in both phases of the theory. This is even 
more remarkable since there is an understanding, at least 
in SU{N) (non-abelian) and in Z^ (abelian) gauge the- 
ories that a perimeter law for both of these order pa- 
rameters requires the existence of massless particles [9J . 
It is not unreasonable to believe that non-compact C/(l) 
gauge theory that we study, will behave in the same way, 
especially given that the large N limit of Zjy gauge the- 
ory does give compact U{1) gauge theory. However, we 
have no massless particles in the theory. We believe that 
we manage to circumvent the previous conclusion since 
we have a statistical long range interaction between the 
anyons. 

To see that the order parameters do indeed exhibit a 
perimeter law, we refer to our previous article [1,. In 
this article it was established, through detailed numeri- 
cal analysis, that the average of the Wilson loop exhibits 



a perimeter law as a function of its parameters L and T. 
It was also found that it exhibits this behaviour on ei- 
ther side of the phase transition, which occurs at around 
/i « 0.152. However the Wilson loop average is calculated 
from the linking number of the fixed loop of dimension L 
and T with the dynamical vortex loops. But the expec- 
tation value of the 't Hooft loop is based on exactly the 
same calculation. For the Wilson loop we compute the 
average value 
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while the 't Hooft loop we compute 
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However for a given rectangular loop, A'^h'l — N'tH] 
hence, the perimeter behaviour that was established for 
the Wilson loop in [T goes over to the 't Hooft loop. 
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